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Abstract. Cosmological Gravitational Waves (GWs) are usually associated with the 
transverse-traceless part of the metric perturbations in the context of the theory of 
cosmological perturbations. These modes are just the usual polarizations '+' and 
' X ' which appear in the general relativity theory. However, in the majority of the 
alternative theories of gravity, GWs can present more than these two polarization 
states. In this context, the Newman-Penrose formalism is particularly suitable for 
evaluating the number of non-null GW modes. In the present work we intend to 
take into account these extra polarization states for cosmological GWs in alternative 
theories of gravity. As an application, we derive the dynamical equations for 
cosmological GWs for two specific theories, namely, a general scalar-tensor theory 
which presents four polarization states and a massive bimetric theory which is in the 
most general case with six polarization states for GWs. However, the mathematical tool 
presented here is quite general, so it can be used to study cosmological perturbations 
in all metric theories of gravity. 



1. Introduction 

The future detection of gravitational waves (GWs) of cosmological origin will strongly 
constrain the possible inflationary scenarios which have been proposed in the last 
decades. Also, GWs will be useful to distinguish between the standard inflationary 
model and the alternative early Universe cosmologies, like the Pre-Big-Bang scenario, 
since the predicted power spectrum of each model can present very different features. 

In the general relativity theory (GRT) the usual procedure in order to evaluate the 
power spectrum of cosmological GWs is to expand small metric perturbations around 
the spatially homogeneous and isotropic Friedmann-Robsertson- Walker (FRW) metric. 
The next step is to identify GWs with the transverse-traceless (TT) part of the metric 
perturbations which do not couple with the perturbations of the perfect-fluid. Thus, 
once it was generated, this radiative gravitational fleld can freely travel through the 
space and reach an observer today. 

The observational effect of GWs is to generate relative tidal accelerations between 
test particles. The Riemann tensor determines these relative accelerations and it is the 
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only locally observable imprint of gravity. To see that, consider a freely falling observer 
at any fiducial point P in the region. Let the observer set up an approximately Lorentz, 
normal coordinate system {a;^} = {t,a;*}, with P as origin. For a particle with spatial 
coordinates at rest, the relative acceleration with respect to P is: 

ai = -RiQjox\ (1) 

where -Rjojo are the so-called 'electric' components of the Riemann tensor due to waves 
or other external gravitational influences. When the linearized theory is considered, the 
Riemann tensor can be split in six algebraically independent components, but for the 
vacuum field equations of GRT they reduce to two, which represent the two polarization 
states of free GWs. These are the 'TT-modes', also called the + and x polarizations. 

Although these are the most studied modes in the GWs physics, when the 
framework of an alternative theory of gravity is considered, the number of non-null 
components of the Riemann tensor can be greater than two and the theory presents 
not only the TT-modes, but also other polarization states can appear. This is a 
direct consequence of the new field equations which can generate other radiative modes. 
Thus, for a generic theory of gravity, GWs can present up to six polarization states 
corresponding to the six independent components of -Rjojo- 

Therefore, we should state that in order to work out the cosmological metric 
perturbations in an alternative theory of gravity, the first step is to find the number 
of independent polarizations of GWs in such a theory, i.e., the number of non-null 
components of the Riemann tensor. This can be done by a very elegant method which 
consists in the evaluation of the non-null Newman-Penrose (NP) quantities Pfl |2] of a 
given theory. These quantities are the irreducible parts of the Riemann tensor written 
in a complex tetrad basis which make them very useful to evaluate the polarizations of 
GWs in an unambiguous way. 

Thus, the present paper intends to consider the general formalism of cosmological 
perturbations in the context of alternative theories of gravity, focusing on the dynamical 
equations of the GW modes. The theory of cosmological perturbations in GRT has been 
largely studied in the literature. Some classical examples are the works by Lifshiftz [3], 
Bardeen [3], Peebles [5], and Mukhanov, Feldman and Brandenberger [6]. For a recent 
review on cosmological dynamics see, e.g.,[71|8]. In the usual approach, the TT-modes 
of GWs are consistently described by the superadiabatic, or parametric, amplification 
mechanism in the GRT p]. Furthermore, it was shown by Barrow and de Garcia Maia 
[Tin E] that the same mechanism applies for modified theories of gravity as scalar- 
tensor theories and the so-called f{R) theories. Although they have analyzed only the 
evolution of the two TT-modes, it is known since the work by Eardley et al. [Tj that 
scalar-tensor theories present in addition at least one scalar GW polarization (more 
general scalar-tensor theories present two scalar GW polarizations) as a consequence 
of the additional degree of freedom included by the scalar field. The relic scalar GW 
production which arises from this kind of theories was recently discussed by Capozziello 
et al. [12], and an upper limit was obtained from the amplitude of scalar perturbations 
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in the Wilkinson Microwave Anisotropy Probe (WMAP) data. But in their analysis, 
Capozziello et al. have considered only the vacuum field equations, this is a limitation 
of their results since a coupling between scalar GWs and the scalar perturbations of the 
cosmological perfect fiuid are expected as will be clear in our derivations. 

In the case of f{R) theories, using the NP formalism, it was shown that a particular 
class of functions f{R) presents two scalar GW modes in addition to the + and x modes, 
thus totalizing four independent polarizations of GWs [13] . However, when the Palatini 
approach is used in the derivation of the field equations, the theory reveals only the usual 
two TT polarizations. It was also found that the scalar longitudinal mode which appears 
in /(-R) theories is a massive mode which is potentially detectable by the future space 
GW interferometer LISA [H]. Again by considering only the vacuum equations, the 
production of the relic GWs of this particular mode was also considered and constraints 
using the WMAP data was established [15]. 

Moreover, the study of extra polarization states of cosmological GWs in the context 
of alternative theories of gravity can reveal new interesting features of these theories 
which do not appear in GRT. A remarkable example is the presence of vector longitudinal 
polarization modes of GWs in some theories. These modes give rise to a non usual Sachs- 
Wolf effect which leaves a vector signature on the CMB polarization [16]. Otherwise, 
vector perturbations in GRT decay too fast and it would not leave any signature on 
CMB polarization. 

Therefore, it is clear that the future detection of GWs, and the corresponding 
determination of the number of polarization modes, are powerful tools to test the 
underlying gravity theory. Thus, the goal of the present paper is to furnish a general 
formalism to find the evolution equations of all the possible polarization modes which 
could appear in a generic theory of gravity. Once the number of independent polarization 
modes are found and the corresponding evolutionary equations could be obtained, one is 
in a position to obtain the power spectrum of each mode, finding the CMB signatures and 
constraining the additional modes. In order to show the application of the formalism, 
and in order to find some new features of the current studied theories, we have chosen 
to obtain the dynamical equations for GWs in the context of two particular theories, 
namely, a general scalar-tensor theory and a bimetric massive theory of gravity. 

First proposed by Brans and Dicke [IT] in the aim of making the theory of gravity 
compatible with the Mach's principle, the scalar-tensor theories are of a great interest 
since, as pointed out by several authors, a coupling between a scalar field and gravity 
seems to be a generic outcome of the low-energy limit of string theories (see, e.g., [18]). 
Another interest in the scalar-tensor models is that the f{R) theories can be written 
as the Einstein equations plus a scalar field, and thus we could in principle extend the 
same formalism applied for the scalar-tensor theories to the f{R) field equations. The 
bimetric massive theory we consider was proposed by Visser [19] in the aim to obtain 
general covariant field equations with massive gravitons. His method was based on the 
introduction of a non-dynamical metric {go)f_iu besides the physical metric (7^^. The 
resulting equations appear as a small modification of the Einstein field equations for 
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which the massive gravitons and the metric {go)fiu are present only in an additional 
energy-momentum tensor. Furthermore, our past studies have shown that the Visser's 
theory is a potential explanation for the current acceleration of the expansion of the 
Universe [201 [21]. 

In deriving the equations for GWs in the two theories we will first review how 
to obtain the number of independent polarization modes for any theory following the 
Eardley et al. approach PJ. In the case of the scalar-tensor models the theory present 
four polarization states in the more general case. Otherwise, the Visser's theory is a 
simple example of how a weak modification of gravity can produce six polarization 
modes. The subsequent analysis show that all the polarization modes, apart from 
the usual -|- and x polarizations, are dynamically "coupled" to the perturbations of 
the cosmological perfect fluid. We argue that this kind of coupling and the existence 
of additional polarization states could furnish distinguishable signatures of alternative 
theories in the power spectrum of the relic GWs. 

The paper is organized as follows: in the section |2] we present an overview of the 
NP formalism starting from the definition of the NP quantities which define the six 
possible polarization states for GWs. Then we find the non-vanishing parameters for 
the GRT, scalar-tensor theories and for the Visser's model. In the section [3| considering 
a generic theory, we find general expressions for the perturbed Einstein tensor and for 
the energy-momentum tensor in the generalized harmonic coordinates. In the section |4] 
we introduce a decomposition scheme which depends on the number of non-vanishing 
polarization modes of GWs which could appear in the various alternative theories. In 
the sections [5] and |6] we apply the formalism of the preceding sections for two particular 
theories, the scalar-tensor theory and the Visse's bimetric model. Finally, we present 
our conclusions and discussions in the section [71 

Throughout the paper we use units such that c = 1 unless otherwise mentioned. 

2. An overview of the Newman-Penrose formalism 
2.1. Tetrads components of tensors and null tetrads 

At every point of the space one can introduce systems of four linearly independent 
vectors el^\ which are known as tetrads. The index in parenthesis is the tetrad index 
which numbers the vectors from one to four. We can define the matrix: 



which is an arbitrary symmetric matrix with negative-definite determinant. Its inverse 
9ic){a), which is defined by: 



9 



(a)(fe) 




(2) 



9ic)ia)9^^^^'^ = S\ 



(3) 



can be used to lower the tetrad indices: 



(4) 
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and to solve (|2| for g^^. 

g,. = giamefe^\ (5) 

One can write any arbitrary vector or tensor as a linear combination of the four 
tetrad vectors: 

rj^^iu- _ mia)ib)- M pi'^p(c)„(d) . .. (p.) 

— J- (c)(d)-6(a)*^{6)S '^'^ ' V^J 

where the quantities T^'^^^''^'" (^c){d)- are the tetrad components of the tensor. They are 
calculated according to: 

ic){d)- — J- -ycr-e^ ^{c)^{d) 5 \' ) 

which is consistent with ^ and Tetrad indices are raised and lowered with (yi 
and 5f(a)(b) respectively. 

The advantages offered in many cases by the use of the tetrad components become 
clear when one examines their transformation properties and when one introduces 
tetrads which are appropriate to the particular problem being investigated. From the 
equation ([T]) one can see that the tetrad components behave like scalars under coordinate 
transformations, i.e., the tetrad indices of tensors do not change under a coordinate 
transformation. Therefore, we have a good way of investigating the algebraic properties 
of tensors in a coordinate-independent fashion by the choice of the tetrads. 

A possible choice is to identify the tetrad vectors with the base vectors of a Cartesian 
coordinate system in the local Minkowski system of the point concerned: 

9{am = e1^f'(f,)9f.u = V{a)(b) = diag(-l, 1, 1, 1). (8) 

The four tetrad vectors, which we shall call (e^, e^, e^^, e^), form an orthonormal 
system of one timelike and three spacelike vectors. 

Another special case is the use of null vectors as tetrad vectors. A particular tetrad, 
known as Newman-Penrose tetrad [2J can be constructed from the orthonormal system 
introduced above, two of the four vectors k and 1 are real null vectors: 

k=^(et + e^), \= ^{et-e,), (9) 

and the other two null vectors m and m are compex conjugates of each other: 

^(e. + .e,), m=i_ 

It is easy to verify that the tetrad vectors obey the relations: 

— kl = mm=l fll 



m = — (e^ + ie^), m = — (e^. - ie^). (10) 



km = km = lm = lm = 0, (12) 



9 {am 



we obtain: 
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/ 



(13) 



A null tetrad basis is especially suitable for discussing null or nearly null waves. 
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2.2. Tetrad components of the Riemann tensor 

The Riemann tensor Rx^y^ can be split in the irreducible parts: the Weyl tensor, the 
traceless Ricci tensor and the curvature scalar (see, e.g., [22]), whose tetrad components 
can be named respectively as \E', $ and A following the notation of [2]. In general, in a 
four dimensional space we have ten \E''s, nine $'s and one A which are all algebraically 
independent. However, when we restrict ourselves to nearly plane waves, we find that 
the differential and algebraic properties of R\^uk reduce the number of independent 
components to six [T]. Thus, following Eardley et al. [1] we shall choose the set 
{\E'2, \I/3, \l/4, $22} to describe, in a given null frame, the six independent components 
of a wave in a generic metric theory. These NP quantities are related to the following 
components of the Riemann tensor in the null tetrad basis described earlier: 





(14) 


^2 — —-^Rlklk, 



^3 = —^Rlklm, 


(15) 


^4 = —Rmm, 


(16) 


$22 = —Rlmlm- 


(17) 



Note that, \E'3 and \E'4 are complex, thus each one represents two independent 
polarizations. One polarization for the real part and one for the imaginary part, thus 
totalizing six components. Three are transverse to the direction of propagation, with two 
representing quadrupolar deformations and one representing a monopolar "breathing" 
deformation. Three modes are longitudinal, with one an axially symmetric stretching 
mode in the propagation direction, and one quadrupolar mode in each of the two 
orthogonal planes containing the propagation direction. The Fig. [T| which was taken 
from [23], shows the displacements induced on a ring of freely falling test particles by 
each one of these modes. GRT predicts only the first two transverse quadrupolar modes 
(a) and (b). 

Other useful expressions are the following relations for the Ricci tensor: 



Rlk — Rlklk, 


(18) 


Rll = '^Rlmlmi 


(19) 


Rim Rlklmi 


(20) 


Rim Rlklmi 


(21) 



and for the curvature scalar: 

R = —2Rik = —2Rikik- (22) 

The overall relative accelerations in a sphere of test particles is described by the 
relation ([T]) and can be expressed in terms of the symmetric "driving-force matrix" 
S whose components are not but the electric components of the Riemann tensor 
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Sij = Riojo, where the latin indices represent spatial coordinates. These components 
can be written as a combination of the NP quantities in the following way: 



/ -|(Re^4 + $22) 

^Im*4 
-2Re^3 



^Im^f4 

|(Re^4-$22) 
2Im^, 



-2Re^3 
2Im^'3 

-6^2 



\ 



(23) 



Gravitational-Wave Polarization 



(a) 


o 










(c) 















e 


(b) 






Figure 1. The six polarization modes of weak, plane, null GW permitted in any 
metric theory of gravity. Also shown is the displacement that each mode induces on a 
sphere of test particles. The wave propagates out of the plane in (a), (6) and (c), and 
it propagates in the plane in (d), (e) and (/). The displacement induced on the sphere 
of test particles corresponds to the following NewmanPenrose quantities: ReV['4(a), 
Im*4(fe), $22 (c), *2(d), Re^-sle), Im^3(/). See ref. [53]. 

Since each NP amplitude is linearly independent, we can expand the components Sij as 
the sum: 

Epw4^ (24) 



1-3 



r=l 



where we have renamed the NP quantities as follows [T]: 

ViX) = ^2, 



(25) 
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P(2) = Re*3, 

P(4) = Re*4, 
P(5) = Im*4, 

P(6) = $22, 



(26) 
(27) 
(28) 
(29) 
(30) 



and E'lf are the components of the "basis polarization matrices" which are given by:. 



/ \ 



E(i) = -6 




1 



E 



(2) 



-2 



/ 





1 



E 



(3) 



/ \ 
1 
1 



E 



(4) 



/l 





\ 

-1 




(31) 



E 



(5) 



/ 1 \ 
1 




E 



(6) 



/ 1 \ 
1 




Finally, analyzing the behavior of the set {\E'2, ^3, ^4, '^'22} under rotations, we 
see that they have the respective helicity values s = {0,±1,±2,0}. These are all the 
possible hehcity values for GWs in a general theory of gravitation. 



2.3. Determining the polarizations of GWs for some specific theories 

The procedure of evaluation of the number of independent polarizations involves 
examining the far-field, linearized, vacuum field equations of a theory, and then finding 
the non-null NP amphtudes. For GW sources sufficiently far from the observer we hope 
that the external solutions of GWs approach that of vacuum and linearized regime, 
in this way the method becomes an unambiguous tool for the determination of the 
propagating physical modes of GWs. In what follows we are going to evaluate the NP 
components of the Riemann tensor for three theories, namely, the GRT, a general class 
of scalar-tensor theories and a bimetric massive theory. 

In the case of the GRT, the field equations can be obtained from the Einsten-Hilbert 
action: 

^ ^ — - - (32) 



—gRd X + Im, 



IGttG 

where Im is the action which describes the matter fields. From the Hamilton principle 
SI — and the Einstein equations read: 



G 



-SttGT. 



(33) 
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For the vacuum T^j, = and (33) reduces to: 

R^.u = 0. (34) 
Therefore, from the relations for the Ricci tensor and for the curvature scalar in 



the NP tetrad basis (18)-(22) we can show that: 

Rlklk = Rlklrh = Rlmlm = 0, (35) 

or equivalently: 

^2 = ^3 = ^22 = 0, (36) 

and since we have no further constraints on the components of the Riemann tensor we 
conclude that: 

^4 ^ 0. (37) 

Hence, as expected, a GW in the GRT presents two polarization states with hehcity 
s = ±2. 

The scalar-tensor theories were first proposed by Brans and Dicke [17] in the aim of 
making the theory of gravity compatible with the Mach's principle. These theories are 
of great interest since a coupling between a scalar field and gravity seems to be a generic 
outcome of the low-energy limit of string theories (see, e.g., [E]). Here, we will consider 
a class of theories characterized by the general action discussed by Wagoner [21] and 
Nordtvedt [25], which we write as found in Maggiore and Nicolis [26] : 

/ = y V^d^'x [-ipR + ip-'uj{ip)VaW'ip - U{ip)] + Im, (38) 

where ip is a scalar field, u^ip) is the coupling parameter and U (ip) can be interpreted as 
a potential associated with ip. The Brans-Dicke action can be obtained for the special 



case when u is constant and U{(p) = 0. The field equations obtained by minimizing (38) 
with respect to the metric and with respect to the scalar field are: 

Gf,u = ~~^/^^~^ - ^9f,u^;a'P''^^-^{^;f,u-gf,uOLp)-Ugf,u, (39) 

[3 + 2tuiip)] Dip-ip^ = 8nT - ^¥^.<,</^'", (40) 
dip dip 

= 0. (41) 

Note that the last condition, which express the conservation of the energy- 
momentum tensor, must be independently imposed since it is not a direct consequence 
of the field equations as in the case of the GRT. The vacuum field equations now read: 

R^lu = -'^^^ip-^.f-u - ^ (^^■,f,u + ^g^^uOip^ + gf,uU{ip), (42) 



and 



R = -^ip-a^'''-3^+AU. (43) 
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Let us first evaluate the polarizations for the case of the Brans-Dicke theory, for 
which the equation (40) becomes simply: 

□yp = 0, (44) 

with the solution: 

<^ = <^o + <^le^''""^ (45) 

where (pQ is a constant obtained from the cosmological boundary conditions, ipi is a 
small amplitude in such a way we can work only to first order in yji, and is the wave 
vector which is null for this particular case. It follows that the curvature scalar is null: 

= ^ Riuk = 0, (46) 

and the Ricci tensor takes the form: 

R^.u = -— e^^'^^'^g^g.. (47) 

Hence if the wave is propagating in the +z direction, the only non-null components 
of the Ricci tensor are R^z-, Rzt and Ru-, which means that Ru is the only non-null 
component in the tetrad basis, leading us to conclude that: 

Rikim = and Rimim 7^ 0, (48) 

and therefore: 

\l/2 = vE^g = 0, \l/4 and $22 7^ 0. (49) 

That is, GWs in the Brans-Dicke theory presents the two s = ±2 polarizations of 
GRT plus a perpendicular breathing polarization mode with helicity s = 0. 

Otherwise, if the potential U{(f) is not null we have a more general scalar-tensor 
theory. Now, dip is not null in general, but obeys the following relation in the linearized 
regime p6] : 

Dip - mlip = 0, (50) 

where: 

2 "PoidU/dip) 

° 3 + 2Lo{ipo) ' ^ ' 

and !p>Q is the value of (p for which the potential minimum is evaluated. The solution of 
(50) is given by (45), but now the wave vector satisfies: 

g„g° = -ml, (52) 

thus the curvature scalar does not vanish but gives: 

i? = -Sm^^e*"""", (53) 

for first order in (pi, and the equation for the Ricci tensor is: 

(54) 



R = -—€ 
^0 



1 2 
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Therefore, from the evaluation of the tetrad components of the Ricci tensor and 
from the preceding relations for the Riemann tensor, we conclude that for a general 
scalar-tensor theory: 

^3 = 0; ^2, ^4 and $22 7^ 0. (55) 

Now, we have a longitudinal scalar polarization mode besides the three modes which 
appear in the Brans-Dicke theory. Note also that the appearance of the longitudinal 
mode is related to the presence of a massive scalar field in the theory. 



The next theory we are going to analyze is an alternative theory of gravity which takes 
into account massive gravitons by the action p!U]: 



d X 



+ m Cmass{9,9o) 



+ 1 



Ml 



(56) 



where m is the graviton mass in natural units and the Lagrangian Cmass is a function 
of the physical metric g and of a non-dynamical prior defined metric g^. The massive 
Lagrangian proposed by Visser is: 

1 



(57) 



1 2^ 

^{9- 9o)pu - 2 [i9o'T"{9 - 9o),.u] |. 



By varying the action (56) with respect to the physical metric we obtain the field 
equations: 

1 



where the massive tensor M^^" reads: 



M^'' = {go'r i9-9oU-M.pi9o'r'i9-9oU {9o'r ■ 



(58) 



(59) 



Among the existing bimetric theories of gravity (Rosen's bimetric theory [27], for 
example), the most common criterion for the choice of the non-dynamical background 
metric is to impose the Riemann-flat condition R^^vk{s,o) = [231 EB]- Thus, the 
simpler choice is a fiat metric which recovers the Minkowski metric by going to cartesian 
coordinates [121 [SO]- By introducing this background metric, Visser has constructed the 
Lagrangian (57) aiming a general covariant description of massive gravitons in such a 



way that the model circumvents the van Dam-Veltman-Zakharov discontinuity [291 EQ] , 
a inconsistency which plagues other massive terms. 

Furthermore, the Visser's model could be an alternative explanation to the current 
acceleration of the expansion of the Universe as indicated by the recession of distant 
Supernova [21]. In this sense, the massive tensor M^i, mimics the dark energy effects in 
large scales while the theory passes all the local tests of gravity. Hence, it is a matter of 
interest to distinguish between the Visser's theory and the Einstein gravity by dynamical 
tests such as GW observations. 
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From (58) we see that the vacuum field equations with massive gravitons read: 



R^"" = ( M^^ - ^r^'^^M ) , (60) 



and for the weak field approximation the right-hand-side of (60) is simply m hni, where 



hf^^ is a metric perturbation around the fiat background metric g^i, = {go)fiu + h^y with 
\hfiv\ 1. Therefore, since there are no further restrictions on R^y, we conclude that 
all its components are non-null and we find that: 

Rlklk, Rlmlm, Rlklm and Rlklfn 7^ 0, (61) 

and then: 

^2, ^3, ^4 and $22 ^ 0. (62) 

Thus, GWs in the Visser's theory present all the six possible polarization states 
showed in the Fig. [l] This result was first obtained by de Paula et al. [31] in a different 
but equivalent approach. Note that in the limit m — )■ 0, the polarization modes \E'2, 
and $22 vanish and we recover the GRT with the only non-null mode \l/4. 



3. Cosmological Perturbations without decomposition 

Let us consider a general theory of gravity for which the field equations can be written 
in the form: 

G^, + F^, = -SnGT^,, (63) 

where besides the Einstein tensor Gf^^, and the energy-momentum tensor for the matter 
fields Tfj^i,, we have a general tensor F^j, = F^i,(5f"^, 7"^, cci", .) which could be a 
function of the physical metric g^^^, of some prior defined metric •j^^, of vector fields 
zUfj,, of scalar fields ip and of derivatives of these quantities. A prior defined metric is in 
general considered in bimetric theories of gravity, for which besides the dynamical metric 
gfj^u, a kind of "absolute" geometry is specified through 7^^. One of the most known 
theory of this kind is the Rosen's theory ^27j for which the second metric takes into 
account the effects of inertial forces. Another example is the bimetric massive theory 
considered by Visser which was introduced in the last section and we will analyze in 
more detail in the section [6l 

Supposing that there is a cosmological solution of such a theory, we can work out the 
perturbations on a cosmological background metric (7^1,. We adopt the metric 5^^^ as the 
Robertson- Walker metric written in Cartesian coordinates with zero spatial curvature 
k = 0. With these considerations the line element reads: 

ds^ = a\ri)rif,ydx^'dx'', (64) 

where the scale factor a{r]) is a function of the conformal time x° = rj. The cosmic time 
t is related to the conformal time by the relation a{ri)dri = dt and the Minkowski metric 
is 77^^ = diag(-l, +1, +1, +1). 
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Perturbing the metric (64) in the general form 
we now have: 



a'^{v)h,j.u{x'^) with \h^^\ < 1, 



(65) 

In this form, the indices of /i^i, are raised and lowered by the metric rj^,^. With the 



line element (65) in the equations (|63j) we can obtain the perturbed field equations: 

6G^, + 5F^, = -Stt [G5T^, + 5GT^,] , (66) 

where, in order to take into account the theories with varying Newtonian "constant", 
we have included the perturbation 5G. 

Hereafter we will use the generalized harmonic coordinates discussed by Bicak and 
Katz [32] ■ For this coordinate system, we have the condition: 

0, (67) 



where ^F^^ is the metric connection perturbation. It is easy to show that, the condition 
(67) is equivalent to: 

V,(5^^^ = 0, (68) 

where 5g^u = Sg^^ — g^^6g/2 with 6g = g^^^Sgap- One can show that 5g^i, = a^h^^ where 
the hau is the trace-reverse perturbation: 



-h. 



(69) 



with h = rj^'^h/j^u and h 

With the metric perturbations given by (65) and with the condition (68) we can 
evaluate the components of the perturbed Einstein tensor 6G'j^ for first order in h^i^. A 
straightforward calculation leads to: 

1 



SG'n 



5GI 



2a2 



1 

2^ 



1 

2^2 L 



in'^ -'H')hi + Aa'Hdihi 



+ 2a-^'H7]'^{dfhl-dJi) 
- + ^% - + in" + H')hl5i 



(70) 



(71) 



(72) 



The prime in the above expressions denotes derivatives with respect to the 
conformal time ?7, and we have defined the Hubble parameter for the conformal time as 
l-i = a' /a. 

Now, it is necessary to evaluate the perturbations of the energy-momentum tensor 
for a perfect fluid: 



T^, = {p + P)U^U, + Pg^ 



(73) 
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where p and P are the energy density and the pressure, and Uv = —U" = (1,0,0,0) 
is the fluid four-velocity. Considering first order perturbations in each one of these 
quantities the perturbed energy-momentum tensor reads: 

T^. ^ f^u = T^u + ST^u, (74) 

where: 

5r/ = (p + P)g^''{U^6Ux + 6U^Ux) + {6p + 6P)U^U'' + 5P(5/ 

-{p + P)6g^''U^Ux. (75) 

And considering the perturbed metric as defined earher we find the components of 
the perturbed quantity ST^": 

6T^ = -6p, 6T^ = ip + P)iW + a-%), 6T^ = 6P61 (76) 

where, for completeness, we have considered V = 6U'' ^ 1 as the spatial part of the 
perturbation of the fluid four- velocity as it appears in the references [HI E3] : 

U^ = U'' + 6U' = {l-^6goo, V'). (77) 

It is easy to see that the four-velocity is approximately the unit timelike vector 
since we assume l^* <^ 1, and terms proportional to and Vh can thus be neglected. 

Then, from the calculation of the first order perturbation of the tensor SFfj^^, and 
using the perturbed Einstein tensor and the perturbed energy-momentum tensor given 
above, one can write the perturbed field equations for a generic theory (see eq. 66). 

As we will see in the following sections, the above equations are written in a 
suitable way to consider extra-polarization states of the GWs in any alternative theory 
of gravity with the field equations (63). They are also useful to find solutions without 
any decomposition as it was presented in the interesting work by Bicak et al. |33]. Note 
that, despite of the difference in notation, a similar set of equations can be found in 
that reference. 



4. Metric perturbations with extra polarization states of GWs 

Instead of considering the minimal decomposition which appears, for example, in the 
reference [0] , let us write the components of Sgfj^^y in a more general form which was first 
presented by Bessada and Miranda [T6] : 

where (p and B are scalar perturbations and Si are the components of a divergenceless 
vector perturbation diS'^ = 0. In our approach, the quantities 0, B and Si are dynamical 
quantities which appear in any theory. On the other hand, the way the perturbation hij 
must be decomposed depends on the number of independent GW polarization modes, 
and hence, it is theory dependent. Thus, the first step is to evaluate the number of 
non-null NP quantities and then we can decompose hij identifying the GW amplitudes 
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in the expansion. In order to introduce the method, let us first consider the most general 
case for which all the NP quantities are non-null. For this particular case, hij can be 
expanded in terms of six components which represent the six metric amplitudes of GWs: 

6 

(79) 



r=\ 



where ejlJ are the polarization tensors. 

In what follows, without loss of generality, we can consider the wave vector of the 
GW oriented in the +z direction. Thus, we can construct the six 6^*2 by combinations 
of the three ortonormal vectors: 

e = (1,0,0) 

m' = (0,1,0) (80) 

= (0,0,1), 
in the following way: 



(81) 
(82) 
(83) 
(84) 
(85) 
(86) 

One can verify that the tensors (81) - (86) are linearly independent and form an 



(1) 

ij 
(2) 

ij 
(3) 

ij 
(4) 

Jj 
-(5) 

Jj 
-(6) 



ortogonal basis. Writing in a matricial form we have: 
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We can see that, apart from a constant, these are just the basis polarization matrices 



(31) which we have described earlier. Now, the expansion (79) can be written in a more 



intuitive form: 

h 



~l~ "Tij ~\~ Xij ~l~ "^iji 
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where we have defined four new tensors following their helicity values s: 

nij = h\f s = 



' - ■ ^ (89) 



And from the very definition of the quantities (89 ), we have the following properties: 
Xi = ^ = 0, d,^ = = 0. (90) 

Their associated helicity values can be evaluated from the behaviour of the 
perturbations under rotations. Notice also that we have chosen to define separately 
the two tensors for the s = modes of GWs since one of them is a longitudinal mode 
(TTjj) and the other is transversal to the direction of propagation {ipij). Furthermore, 
we have seen in the section [2] that these two modes can appear for some theories (in 
a general scalar-tensor theory, for example), but other theories have a structure such 
that only one scalar mode appears (this is the case of the Brans-Dicke theory). This is 
a particular feature of the scalar modes, since it is not possible that a certain theory 
presents only one of the two tensor modes and not the other. In fact, there is no 
theory for which the + and x polarizations appear separated for vacuum GWs as it 
was evidenced by the construction of the NP amplitude (see section |2|. The same 
happens with the two vector modes which generate the same NP quantity \I'3. 

It is instructive to write the NP quantities in terms of these metric perturbations 
considering an observer today located at the local Minkowskian space-time. Thus, 
considering the Cartesian tetrad basis we have: 

^2 = -^(tT^^.OO - '^T^zz^zz + 4>,zz) + ^(hfio - h^^,), (91) 
^3 = -^^[{'rxz,zz - Tcczfio) + «(Tjy2,00 - Ty2,zz)], (92) 

^4 = ^[XyyfiO - XxxfiO + 2iXxs/,oo], (93) 
*22 = -^{'IpxxfiO + i'yyfiO " ^,Oo), (94) 



where: 



/i = -0 + z^i + 7ri. (95) 



Now, since the polarization modes of GWs are linearly independent, we can 
summarize the procedure to calculate the perturbed field equations for a given general 



theory of the form (63) in the following way: 



Scalar metric perturbations 
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The general scalar metric perturbations are 




(96) 



If the theory has ^2 ^ ^ and $22 7^ 0, it means that an observer today can measure 
the two scalar GW modes (s = 0). Hence, in order to describe the evolution of these 
modes we write the term n't} in the form: 



If the theory has \E'2 = and $22 7^ 0, no longitudinal scalar GW amplitude can 
be measured. Then vrjj is suppressed, and in order to do not change the number of 
degrees of freedom of the metric perturbations, a scalar longitudinal component, say 
Z), must be added. This new scalar component represents a djTiamical perturbation of 
the cosmic gravitational potential rather than a radiative GW field. Consequently, an 
observer today measures only one scalar GW mode whose evolution can be described 

is) 

by ipij, and hi, takes the form: 



where we have introduced D using the spatial partial derivatives in order to represent 
the longitudinal behaviour of this quantity. This is the same procedure of introducing 
longitudinal scalars in the minimal decomposition [6]. For those theories which can 
be found in this case (e.g., the Brans-Dicke theory), the only non-null scalar degree of 
freedom (in the local Minkowskian frame) is ipij. In this sense, it is always possible 
to make D = by gauge transformations in the reference frame of the local far 
field observer, although in general D does not cancel out when considering the whole 
cosmological evolution of the perturbations. 

On the other hand, if the theory has '^2 = $22 = 0, the perturbations must be 
decomposed in the usual way, since there are no scalar GWs which could be observed 
today, that is: 



where we have added the new scalar quantity E to describe the evolution of the 
perturbations of the gravitational potential. The most important theory which can be 
found in this case is the Einstein theory for which we have only the two "pure tensor" 
degrees of freedom for vacuum GWs. In Minkowski coordinates, it is always possible to 
cancel out all the scalar components by gauge transformations, but in the presence of 
the cosmic fluid they are important to describe the evolution of the perturbations. 

• Vector metric perturbations 




hf> = d,d,D + ifj, 



(98) 



(99) 
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The general vector metric perturbations for any theory are: 

59^;}=a\v){ I %\ (100) 



where the vector Si and hlj satisfies: 



i{v) 



diS' = = (101) 

If the NP quantity is non-null (the Visser's theory is an example), an observer 
today can measure the two vector modes (s = ±1) of GWs, and the metric perturbations 
/ij-J'' can be identified to the corresponding GW amplitude: 

^ = n,, (102) 

with = 0. 

On the other hand, if \I'3 = 0, there are no vector GWs today and we have the 
usual representation in terms of the vector quantity Qi\ 

/iS;) = 9,Q, + a,Q„ (103) 



where, from (101) Qi is divergenceless diQ^ = and again, it is easy to verify that 
the number of the degrees of freedom of the metric perturbations does not change. 
The quantity Qi is a vector perturbation of the gravitational potential for which the 
dynamical equations gives, in general, a decaying mode in the context of the GRT 
(see, e.g. [HE])- As in the case of scalar perturbations we have discussed earlier, it 
is always possible to cancel Qi by gauge transformations in the local Minkowski frame 
when \E'3 = 0. 



• Tensor metric perturbations 

Finally, the tensor metric perturbations are constructed using a symmetric tensor 
Xij which satisfies the constraints: 

= d,xi = 0. (104) 

Thus, the tensor component which corresponds to GWs with s = ±2 is written in 
the usual way in a theory- independent form: 

Counting the number of independent components we have used to construct Sg^^, 
and the number of constraints, we can see that we have four functions for scalar 
perturbations, four functions for vector perturbations, and two functions for tensor 
perturbations. Thus, as expected, we have ten independent components of Sg^^- 
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In the context of the Einstein theory we have already found that \I/4 is the only 
non-null NP quantity, thus it is direct to verify that the metric perturbations have the 
usual minimal decomposition: 

og,u a y^^^ ^^^^ ^ ^^^^^ ^ ^^^^ ^ ^^^^ ^ J , uuDj 

and GWs are described only by the quantity Xij whose evolution equation can be derived 
from SG^f^l = —SttGT^u to obtain the very known result: 

xf + 2nxf - v'xi = 0. (107) 

In the following sections we will exemplify the above decomposition scheme for the 
two other theories which we have treated in the section |2| namely, a general scalar-tensor 
theory and the Visser's bimetric theory with massive gravitons. For each theory, after 
identifying the form of the metric perturbations we will find the dynamical equations 



for GWs in the coordinate system defined by (67). 
5. GWs in scalar-tensor theories 

In this section we will consider perturbations of the general class of scalar-tensor theories 



introduced in the section^ With a glance at the field equations (39) we identify it with 



the general form (63) where G{(f) = ^ and the generic function F^^^, takes the form: 

Ff,u = ^^^^ (^V;i,V;u - ^g^,uV■,a'^'"^ + ^(v';/^^ - Qf^uOLp) + g^uU{ip). (108) 
Thus, disturbing the Newtonian "constant" : 

(109) 

and evaluating the components of the perturbation we can find the perturbed field 



equations (66) which for this case reads: 

5G,, = -— (sT,, - ^-^T,}j + 5F,,. (110) 

As we have already shown in the section |2| the evaluation of the NP parameters for 
the most general scalar-tensor theory lead us to conclude that an observer today would 
measure: 

^2 7^0, ^3 = 0, ^4 7^0 and $22 7^0, (111) 

remembering that \l/2 = if the potential associated to the scalar field is null U{(p) = 0. 
Thus, following the procedure of the last section we can write the perturbations for the 
general case which reads: 

y Si + OiB TTij + Xij + Wij + OiQj + OjQi j 

Now, we can introduce the perturbations given above in the perturbed Einstein 
tensor and in the perturbed energy-momentum tensor calculated without decomposition 
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in the section [3| and we can also evaluate the perturbed quantity Thus, after a 

straightforward calculation we find the dynamical equations which describe scalar and 
tensorial GWs in the context of the scalar-tensor theories: 



Scalar 



b" + 2-H0' - {m^ - m') 



CO (if 



2\!f 



+ 



-Ian' - - 2aHW^B = 167raV"^f - 5p + ^p) - 2Ai, (113) 

2 \ if J 



{d'B)" + mid'B)' - V\d'B) + iU' - 'W)d'B + 



Voj' 



d'B 



2Ua-\2d'(j) + dk-K'") = IQnGa'ip + P){Vf + d'B) + 2A 



25 



;ii4) 



+{H' -n^ 



+ 



1 

+ -w— 

4 v9 



2A.-^' 



(115) 



Tensor 



Xi 



2U 



2, J 



0. 



;ii6) 



In the above equations we defined = nl + ip-^, and V^^ is the component of 
which is parallel to the direction of propagation. The perturbed quantities Ai, and 
A^i^ take into account the perturbation of the scalar field ip and its derivatives. They 
are defined as follows: 



Ai 



A,' 



5uj + m — — + a^ 

if if 



uj—d 



(117) 



118) 



As/' 



+ 



;ii9) 



The equation for tensor perturbations in scalar-tensor theories (116) was studied 



in the reference [TOj. It represents the evolution of free GWs with helicity s = ±2 in 
such theories. These modes do not couple with the perturbations of the perfect fluid 
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and are generated quantum-mechanically due to vacuum perturbations in the very early 
Universe. They are amphfied by the expansion of the Universe due the process known 
as superadiabatic amphfication [9]. In this case, the scalar field (p contributes for the 
cosmological potential which generates the amplification. 

Regarding scalar perturbations, the driven equations for the NP modes \E'2 and $22 
are given by the equations (113), (114) and (115). These set of equations, considerably 
more complicated than the tensorial case, show some new physical features when 
compared with the usual metric decomposition. The most important difference is that 
these equations represent the evolution of radiative fields coupled to the evolution of 
the generalized Newtonian gravitational potential. The radiative fields are represented 
by the quantities and ipij while the generalized Newtonian potential is described by 
the dynamical perturbations (j) and B or some particular combination of them. Thus, in 
some sense, we can say that if the scalar-tensor theory is the "correct" theory we would 
have a new kind of cosmological GW background. On the contrary to the tensorial case, 
this new background is coupled to the matter perturbations of the perfect fiuid. This 
is expected since the usual scalar perturbations are also coupled to the fiuid dynamics. 
But, in the usual sense, GWs would be coupled to the matter only if a component of 
anisotropic stress would be present. 

A direct consequence of such a coupling is that the whole evolution of the density 
perturbations of the cosmic fiuid would depend not only on the evolution of (f) and B 
but also on the evolution of the amplitudes of the scalar GWs vTjj and ipij. Thus, from 
(113), for example, we are lead to conclude that: 

^ = ^(0,i?,aM, (120) 

where we have also included the dependence on the perturbations of the scalar field 
5ip. Notice that, in fact, the dependence of 5p/p with the GWs amplitudes appears as 
the dependence on the trace of the overall contribution of the scalar GW amplitudes 
Ci = K + i^i- Therefore, a complete understanding of the evolution of the GWs with 
helicity s = requires the knowledge of the evolution of the scalar perturbations and, 
similarly, these GW modes affect the evolution of the density perturbations. 

Another issue of particular interest is how the presence of the scalar GWs would 
affect the angular pattern of the CMB. Since the GW amphtudes iiij and iptj now enter 
the geodesic equation for photons, it is expected that they leave a signature on the 
CMB due the so-called Sachs- Wolfe effect, which can be understood as the shift of 
photon frequency along the line of sight. The small fluctuations in the CMB may be 
conveniently described by perturbations of the temperature parameter T in the Planck 
distribution /. The computation of the contribution of the scalar GWs to the angular 
temperature inhomogeneities ST/T of the CMB is out of the scope of the present paper 
and a rigorous treatment will appear elsewhere. 
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6. GW modes in a bimetric theory of gravity 



Now, let us turn our attention to the massive bimetric theory first considered by Visser 



[19j which we have introduced in the section |2| Comparing the field equations (58) with 
the generic form (63) we identify F^'^ = —m?M^^. 



Our explicit calculations of the section [2} and the previous result by de Paula et 
al. [31], have led us to the conclusion that all the NP quantities are non-null for the 
Visser's model: 

^2^0, ^3 7^0, ^4^0 and $22^0. (121) 

Thus, following the procedure of the section |4| the perturbations now should be 
written in the form: 



Si + diB -Kij + Tij + Xij + ^ij 



(122) 



With (122) in the perturbed field equations (66), calculating the components of 



5Fau and with the help of the perturbed components of the Einstein tensor and of the 



energy-momentum tensor calculated in the section (|3j), we obtain the perturbed field 
equations for Visser's theory for each group of perturbations: 



Scalar 

0" + 2H(t)' - V^(f) - {m^ - m' - m^a^)(f) 

-[{W - n^) - m^a^{a^ - l)]Ci - AaUV^B = -lQnGa^6p, (123) 

{d'By + AU{d'B)' - V\d'B) + ]^[2{W - H^) + m^a\3 - a^)]d'B 
-2Ha-\2d'(j) + dkir'^) = IQnGa^p + P)(W + d'B), (124) 



e + - V'^ - {W + H')i5l + \m\\a' + 1)^^ 

+ma{did^B) - ^[2{H' - U^) - m^a\a^ - = IGTrGa^SPSj. (125) 



Vector 

S'" + mS" - V^S' + ^[2{H' - H^) - m\\a^ - 3)]S' 

+2Ua-^dkT''' = 16nGa\p + P)(\/_[ + S'), (126) 

rf + 2-HTf ' - VV/ + ^m'a\a^ + l)r/ + Anari'^d^^iSk) = 0. (127) 
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Tensor 



J' 



2, J 



1 



0. 



(128) 



(129) 



Similarly to the tensorial equations for the GRT and for the scalar-tensor theory, 
this last equation describe the evolution of free GWs related to the NP mode '^i. 
But now it is the new term which contains m? that contributes to the parametric 
amplification of GWs. 

The equations for the scalar modes (s = 0) are now given by the set (123), (124) 
and (125). The argumentation for the scalar modes TTjj and ipij is similar to the case 
of the scalar-tensor theories, except for the absence of the scalar field perturbation 5ip. 
Again, the three equations must be solved simultaneously in order to find the evolution 
of the GWs with helicity s = and to find the evolution of the density perturbation 
which is related to the metric perturbations through the equation (123): 

P P 

Conversely, we have now GWs with helicity s = ±1 which correspond to the mode 
^^3. The equation which describes the evolution of the GW amplitudes for this mode is 
the equation (127). Note the similarity of this equation to the equation (128), except 
for the presence of the term containing S** in the equation (127). The presence of this 
term makes the vector GW modes coupled with the vector perturbations since 5* is 
coupled with the fluid vector perturbations through equation (126). Furthermore, the 
perpendicular part of the vector perturbation (which is a pure vector) is a function of 
the quantities and r*-' : 

VI = VI{S\t'^). (130) 

Regarding the CMB anisotropy, the presence of the longitudinal vector modes 
of GWs do not yield the well know version of the Sachs- Wolfe effect which appears 
in the GRT or in the scalar-tensor theory [M]. Theories which present vector GWs 
(\E'3 7^ 0) give rise to a nontrivial Sachs- Wolfe effect which leaves a vector signature of 
the quadrupolar form Y2^±i on the CMB polarization (see detailed discussion in [16]). 



7. Final Remarks 



In the present work we have studied the evolution equations of cosmological GWs in 
alternative theories of gravity. Since the most part of the alternative theories present 
more than the two usual + and x polarizations of GWs, we have addressed the problem 
of how one could take into account the new polarization states in the cosmological metric 
perturbations. 

First of all, we have presented an overview of the NP formalism since it is 
particularly suitable for evaluating the number of non-null GW modes of any theory. 
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Then, we have proposed that the construction of the metric perturbations for a given 
theory should depend on the number of non-null NP parameters of the theory. 

The formalism developed here is quite general and can be applied for a wide range 
of alternative theories of gravity. In order to show that, we have evaluated the evolution 
equations for GWs for two different theories: a class of scalar-tensor theories, for which 
Brans-Dicke theory is a particular case, and a massive bimetric theory. In the first case, 
the theory presents two scalar modes (s = 0) and two tensor modes (s = ±2) of GWs. 
In the case of the bimetric theory, GWs have in addition the two vector modes (s ± 1), 
totalizing six polarization states of GWs, i.e., the most general case in the context of a 
four-dimensional theory of gravity. 

A qualitative analysis of the governing equations of the scalar perturbations have 
shown that the evolution of the density perturbations of the cosmological fluid depends 
not only on the generalized Newtonian potential but also on the amplitudes of the scalar 
GWs. Another direct consequence of the presence of TTj^ and ijjij is the possible signatures 
that this quantities would leave in the angular pattern of the CMB. Such a signature 
might impose strong limits in the amplitudes of the scalar GWs by the analysis of the 
CMB data. Moreover, a remarkable effect on the CMB which have already been studied 
in the literature [H] is a non-usual Sachs- Wolfe effect which appears due the presence 
of the longitudinal vector GW modes on the geodesic equation for photons. Such vector 
GW modes are present, for example, in the massive bimetric theory analyzed here. 

The detection of GWs is a particularly challenging issue and it may be the final 
answer to the "correct" theory of gravity. If GWs present non-tensorial polarization 
modes as discussed in the present paper, we will have a stochastic cosmological 
background of GWs which is a mixture of all the polarization modes. If, in analyzing 
such a background, scalar and/or vector GWs could be found, the result would be 
disastrous for the Einstein theory. 

The evaluation of the response function of the non-tensor polarization modes for 
interferometric GW detectors was carried out in the references [35] and [36] . Particularly, 
Nishizawa et al. [SS] have found that more than three detectors can separate the mixture 
of polarization modes in the detector outputs. But they have considered only separation 
between the three groups: scalar, vector and tensor modes of GWs. Furthermore, they 
have found that, statistically, the GW detectors have almost the same sensitivity to 
each polarization mode of the stochastic background of GWs. In the work by Corda 
[37] . the dectabihty of a particular polarization was discussed, namely, the longitudinal 
scalar component. It was also shown that the angular dependence of such a mode could, 
in principle, allows discriminating this polarization with respect to that of CRT. 

A positive detection for certain modes and a negative detection for others may 
exclude a particular theory or, at least, establish strong constraints to the alternative 
theories of gravity. But it is important to emphasize that, the confirmation by the 
observation of the number of non-null GW modes is not enough to determine the 
"correct" theory, since a number of theories can have the same number of non-null 
modes. It is also necessary to evaluate the spectrum of GWs for each mode and for each 
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theory. The evaluation of the spectrum can only be done following the cosmological 
evolution of the GWs. In this sense, the equations given in this paper are the first 
step for such a computation. Finally, the analysis of the response of a given mode and 
the evaluation of the spectrum can bring crucial answers for the comprehension of the 
gravity in cosmological scales. 
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